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Eddington-inspired Born-Infeld gravity in three spacetime dimensions is studied. We obtain 
analytical and numerical solutions for the scale factor in a three dimensional Friedman-Roberston- 
Walker line element with (i)pressureless dust and (ii) a perfect fluid with p — ^, as matter sources. 
Further, we also find a family of static, circularly symmetric analytical line elements with isotropic 
pressures. For FRW cosmology, we note that if the new parameter k which arises in the theory, is 
positive (k > 0), the solutions are singular (except for the open universe, with a specific condition), 
whereas if «: < 0, they represent non-singular spacetimes. On the other hand, in the circularly 
symmetric static case, our solutions are non-singular for both «; > and k < 0. 
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^ I. INTRODUCTION 

<D , 

, Theories of gravity different from General Relativity (GR) have been actively pursued by many, for a variety of 
' reasons. One such reason relates to the possibility of avoiding the singularity problem in GR involving the occurrence 
[ of a big bang in cosmology or black holes in astrophysics. In a classical metric theory of gravity one is aware that these 
singularities are inevitable, as proved through the Hawking-Penrose theorems under very general and physical 
assumptions. However, it is quite possible that in an alternate theory one might obtain non-singular geometries (for 
qh< example non-singular Friedman-Robertson- Walker type cosmologies) as solutions~a feature which does not exist in 
I ' the FRW cosmology based on GR. It may be noted that the removal/resolution of a singularity is also expected to 
be a basic feature of a quantum theory of gravity. 

In this article, we look into one such alternate theory. Historically, its origin goes back to Eddington who showed us 
how de Sitter gravity could be obtained using an Einstein-Hilbert like action where ^J—gR (Ricci scalar) is replaced by 
—det{Rij) [Ij. Eddington's formulation had the advantage of choosing the connection as the basic variable, instead 
of the metric -therefore, it is essentially an affine formulation. However, coupling of matter remained a problem in 
' this formulation. 

\^ \ We arc also aware, in electrodynamics, of Born-Infeld theory Q, which was introduced in order to get rid of the 
' infinity in the field at the location of the charge/current. A gravity theory in the metric formulation inspired by 
' Born-Infeld electrodynamics was suggested by Deser and Gibbons Q ■ Later VoUick Q worked on the various aspects 
■ of the Dcser-Gibbons proposal in a Palatini formulation and also introduced a non-trivial way of coupling matter in 
CO ' such theories. More recently, Banados and Ferreira Q have come up with a formulation wherein the matter coupling 
', is different and simpler from that introduced in VoUick's work [5| . We will focus in the theory proposed in @ and call 
^ ■ it as Eddington-inspired Born-Infeld (EiBI) gravity, for obvious reasons. Note that the EiBI theory has the feature 
that it reduces to GR, in vacuum. 

It may be noted that the theory wc consider falls within the class of bimetric theories of gravity (also called bi-gravity). 
The current bimetric theories have their origin in the seminal work of Isham, Salam and Strathdee Numerous 
papers on varied aspects of such bimetric theories have appeared in the last few years. The central feature here is 
the existence of a physical metric which couples to matter and another auxiliary metric which is not used for matter 
couplings. One needs to solve for both metrics through the field equations. 

Let us now briefly recall Eddington-inspired Born-Infeld gravity. Since we deal with three spacetime dimensions in 
this article, we prefer to write down the action and ensuing field equations in three dimensional spacetime. The action 
for the theory developed in @, is given as: 



5-H 
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X 
S3 



\gij + nRij I - 



5M(ff,f) (1) 



where A = Variation w.r.t F, done using the auxiliary metric = gij + nRij^q) gives 

q,j = grj + uRtjiq) (2) 
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Variation w.r.t gij gives 

^f^q'^ = ^\f^9'^ - Ky/^T'^ (3) 

In order to obtain solutions, we need to assume a gij and a qij with unknown functions, as well as a matter stress 
energy (T'^). Thereafter, we write down the field equations and obtain solutions using some additional assumptions 
about the metric functions and the stress energy. 

Quite some work on various fronts have been carried out on various aspects of this theory in the last couple of years. 
Astrophysical aspects have been discussed in the references in Q while cosmology in those cited in Other topics 
such as a domain wall brane has been analysed in [lo| . More recently, generic features of paradigms on matter-gravity 
couplings have been discussed in However, in [l2| a major problem related to surface singularities has been 

noticed which has put the theory on shaky ground insofar as stellar physics is concerned. 

Our work here is reasonably modest. In the two subsequent sections, we discuss cosmological and circularly symmetric 
solutions in three spacetime dimensions, successively. In the final section, we briefiy summarize and conclude. Most 
of our solutions are analytical and simple. They also maintain some of the generic features noted in the original work 
of Banados and Ferreira . 



II. COSMOLOGY 

Let us assume a homogeneous and isotropic Friedmann-Robertson- Walker (FRW) line element in 2+1 dimensions, 
given as: 



1 - fcr2 



r^de^ 



(4) 



where k = +1,0,-1 for closed, flat and open universe respectively. The energy- momentum tensor is taken to be that 
of a fluid with T'-' = (P + p)u'^u^ + Pg^'-' . The conservation of energy-momentum leads to 



p=-2-{P + p) 



(5) 



which implies p (x -i^ (for P — 0) and p oc (for P — Further, we assume the auxiliary line element to be of the 
form 



ds^ = -U(t)dt^ + a^{t)V{t) 



dr^ 



1 — kr'^ 



r^de^ 



(6) 



Using the auxiliary metric (g'-' ) , physical metric {g'^ ) and stress-energy tensor (T^^ ) in the field equation obtained by 
varying w.r.t. gij, we get two equations, given by. 



U 



V 



D 



1 + KpT 

D 



(7a) 



(7b) 



1 - kPt 

where, D — {1 + kpt){1 — uPtY , Pt — P + ^, and Pt = P — A. We define two quantities G2{p, A) and F2{p, A), given 
as. 



2;7\ 2k U 



G.(p,A).-(l + C/-- 



(8a) 



F2ip,A) = 1 



k(Pt + Pt)(1 — w — kPt — Kwpx) 
(1 + kpt)(1 - kPt) 



where, in Eq. (j8b[) . we have also assumed P ~ wp. Combining Eq. ([2]) and the results from Eqs. 
Friedmann equation given by. 



H' = 



G2 

2Fi 



(8b) 

|), we obtain the 
(9) 
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Let us now examine two special cases: (i) pressureless dust {w = 0) filled flat universe (fc = 0) with the cosmological 
constant, A = and (ii) radiation dominated {w = i) flat universe with A = 0. 

For the case («), using the equation of state (Eq. ([5])) and the Friedmann equation (Eq. ([9])), we find 



4p(l + p) 



(10) 



where, p — np. Note that from the conservation law pa^ is a constant. Using this in Eq. ()10|) . we obtain the solution 

for the scale factor a{t) for k > as. 



i{t) 



(11) 



where the constant Ci = \p\ > 0. The solution clearly demonstrates that for k > the EiBI theory cannot avoid 
the initial singularity {a(t) can become zero at a finite t and hence we have infinite curvature). However, for k < 0, 
the Eq. ([TU| becomes. 



4H(i-H) 



(12) 



From Eq. (1121) . we note that there exists a maximum density (ps) or, equivalently, a minimum value for the scale 
factor (as). The solution for the scale factor a{t) for k < is given by. 



ait) 



(13) 



where C2 is a constant. It is easy to see that the scale factor is never zero and thus, there is no curvature singularity. 
Both these solutions arc plotted in the top row of Fig. 1. 

For case (M),(i.e radiation dominated universe), the Friedmann equation becomes: 



2p-l 



1 



(1 + P) 



Kk 



(2-p) 



(1 + P)(2-P)^ 
(4-p + 4p2)2 



(14) 



For the flat universe {k — 0), the last term in the square-bracket of the right hand side of the Eq. ([Til) docs not 
contribute and the equation becomes. 



9 2 

-k 



2p- 1 



l-§) (1 + p) 



(l + p)(2-p)^ 
(4- p + 4p2)2 



(15) 



In this case, for k > 0, > for an arbitrary, but physically justifiable value of p (i.e for p > 0). Thus, here also, a 
curvature singularity appears in the solution. Using the equation of state :pa^ = C3, where C3 is a constant, we can 
solve numerically the Eq. (fTSj) for the time evolution of the scale-factor a{t), which is shown in Fig. [T] However, for 
K < 0, the Eq. (fTS]) is rewritten as. 



2 



l + 2|p| 



2 



(|P|-1) 



(l-|p|)(2+|p|)^ 

(4+|p|+4|p|2)2 



(16) 



The presence of the factor (1 — \p\) in the Eq. ([TS]) leads to being negative when |p| > 1 . Therefore, in this case, 

for K < 0, there exists a maximum density or, equivalently, a minimum, non-zero value for the scale factor. 

Let us now define two dimensionless variables fl = (0 < < l),z = Using these we recast the equations as: 



\KH^{n)\ 



1 + 2Q 



(l-17)(2 + f7)2 
(4 + + 4172)2 



(17a) 
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FIG. 1: Plot of the time evolution of the scale-factor a{t) for flat universe {k = 0).The parameters are specified on each frame, 
as is the minimum scale factor for non-singular universe. Ci and C3 are constants as defined in the text. 



Numerically solving the above equation, we plot and note the time evolution of the scale factor using the Eqs. (|17p . 
In the bottom row of Fig. [1] we note a non-singular scale factor for k < -a result showing the existence of a bounce, 
similar to that obtained in 3 -I- 1 dimensions. For k > 0, the solution appears to be singular. The above solutions 
(especially, the ones for the P = case) are instructive because they arc, as far as we know, the only known analytical 
solutions in EiBI cosmology. 

Introducing curvature in the spatial slices (i.e. k = +1,-1) docs not yield anything drastically new in the dust 
(P = 0) case primarily because of the fact that p oc ^.Tlie solutions for k > and k < are given by. 



a(.).J(l-A)c.(^-l),.>0 (18a) 



a(t)^^(l-A)c.(^ + l),.<0 (18b) 

where po is a constant {pa^ — po for P ~ 0). For a. k — +1 solution there is a lower bound on po (po > 1) whereas, 
for k ~ —1, Po is an arbitrary positive, real constant. When fc = 0, we recover the earlier results (Eg . [TT] and Eg. [T3l 
). It is easy to see that there is nothing new in the P = 0, A; ^ solutions. 

However, in the radiation dominated (P = |) case we do get some interesting results, though the main conclusion 
regarding the appearance of a singularity or otherwise, is almost the same. The introduction of spatial curvature 
results in an additional term as shown in the square-bracket of the R.H.S. of Eq. (|14p and this leads to all the 
differences. For a closed universe (fc = +1), instead of a bounce we get an oscillation of the universe for k < 0. For 
K > there is an additional feature implying a maximum value of the scale-factor, along with the singularity. These 
are shown in Figl2] In an open universe (fc = —1), wc do not see any characteristic novelties for k < 0, but for 
K > along with the singularity, we also get, under certain circumstances, a non-singular loitering phase of the early 
universe. This is shown through the plots in Fig. [31 
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FIG. 2: Plot of and the time evolution of the scale-factor a{t) for closed universe (fc — +1). pi is a constant, (pa'^ — pi). 

III. CIRCULARLY SYMMETRIC, STATIC SOLUTIONS 

Let us now turn to a completely different class of line elements -i.e. those which are circularly symmetric and 
static. We consider a simple ansatz for the physical line element gij, 



(19) 



where r{l) is a non-negative function and / extends from minus infinity to plus infinity. r(/) represents the radius of 

a circle at each value of I. The energy-momentum tensor is assumed as, r'^=diag. ^p. Pi, fl^^ Energy-momentum 

conservation, VjT*-' = 0, leads to the requirement. Pi = ^2 = P- Let us further assume the auxiliary line element to 
be of the form 

dsl = -dt^ + u^{l)df + v^{iy{l)d9^ (20) 
where u{l) and v{l) are non-negative functions of I. The field equation obtained from gij variation yields, 

u2 - 1 



P = A, 



A 



The other field equation obtained from T variation yields a single equation given as, 

2' 



u — 1 = —K 



u 
u 



(21) 



(22) 



Trivial solutions can be found by assuming u = 1 which gives r{l) ^ I (flat spacetime) and u = /3 1 and a constant) 
which yields de Sitter spacetime in three dimensions. 

Let us now rewrite the equation using the new variables U and S with v?' ~ U,r'^ ~ S in Eq. p2|) . The above 
equation then becomes 



U -I 





fS' 


U'Y 


S' 












+ 45 







(23) 
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FIG. 3: Plot of and a{t) for open universe (fc = — 1) showing the loitering phase at its early stage, pi is a constant, (pa^ = pi 



To solve this equation we need to assume a relationship between U and S. Let us take 



s ^ u 



D 



(24) 



where Z? is a constant. Choosing D appropriately, we can transform the second order differential equation (Eq. 
into an integrablc ODE and solve it. Though different choices of D may give different mathematical solutions, only 
those are physically important which result in physically acceptable energy-density distributions (pt)- One such class 
of solutions is obtained by choosing D as a step function given by. 



where /i > 0, 7^ 1, and v ^ 0,1. 
For K > 0, we get the solution as. 



2/i 
2v 



D = [6{l) - 6{-l)] for K > 



[Oil) - e{-l)] for K < 



(25) 



S — Sf) 



exp 



2 / 2|i| 



(26a) 



U 



1 — /i-^ exp 



2imi-/^") 



(26b) 



where, we have assumed S{0) — Sq and U(0) — 1. The Ricci scalar (i? for the physical metric) for such a solution is 
as follows: 
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We note that such solutions with k > are non-singular. 
When K < 0, wc obtain the solution as, 



S 



exp 



2fUI \ , 2 

— — ' v'' exp 



2|i| 



(28a) 



U = 



I + ly"^ exp 

For this solution, the Ricci scalar turns out to be, 

2''KI ^ I 1 ( 2|/| 



2|i|(l-H^^) 



(28b) 



i? 



2z.^ 



exp 



IJ^exp 



vy'\K.\ 



exp 



2i/|i| 



exp 



2\l\ 



i/-^|k| 



exp ' — exp 



2\l\ 



uW\k.\ 



exp [ —j= ] + i^^ exp 



2\l\ 



(29) 



For = — 1, Eq. (|28ap becomes. 



S ^ So cosh 



2/ 



(30) 



The corresponding Ricci scalar, from Eg. turns out to be 

2 + 2sech^ I 



R 



which is non-singular. The Kretschmann scalar is 



■ijkl 



4 




2/ Y 




2 - tanh ^ 


W 





(31) 



(32) 



and it is non-singular as well. Therefore, assuming 5*0 — Tq, the physical line element, for this case (k < and 
ly ~ —I), can be written as. 



ds^ = -dt^ + df + rl cosh 



21 



de^ 



One can also rewrite it as. 



where, 



ds' = -dt' 



dr^ 



b(r) 



r^d9^ 



1 /r^ 
bir) =r + —[^-r'' 
K V r 



(33) 



(34) 



(35) 



and r = rQi cosh 



Thus, r > tq and a metric singularity occurs in the transformed metric function only at 



r ~ rQ. We plot the transformed metric function 1/(1 — ^^) and have shown it in Fig. Note that this solution is 

not asymptotically flat (i.e. does not tend to zero as r — > oo). However, there is a minimum radius r = tq and 
the geometry is symmetric as I — > — L Also b{r = vq) ~ tq. Thus, the features are similar to that of a Lorentzian 
wormhole |l3| though the geometry is not asymptotically flat. 
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FIG. 4: Plot of ^ 



b(r) 



In Fig. Owe have shown and compared the g and q metric functions, the Ricci scalar (i?), and the energy-density 
(pt) for two different solutions, one for k > 0, = 2.0 and the other for k < 0, v ~ —1.0. Note that all the plots are 
symmetric as / — > —I. In Fig. [S]we have plotted the Ricci scalar of the g metric and the energy density for different 
solutions. 



In Fig. 6(a) we have shown the solutions for k > with two different values of the parameter fi, one for ^ = 0.5 and 
the other for /i = 3.0. In both solutions, we note that the Ricci scalar is symmetric and tends to a negative constant 
value for large positive or negative /-it is asymptotically AdS, but near the origin it is positive. Interestingly, for 
/i — > 0, or /i 3> 1, it tends to zero (asymptotically flat spacetime). The energy density is also symmetric as Z — > — / for 
both fi = 0.5, 3.0. In fact, by suitably adjusting A one can maintain a positive energy density everywhere for all /. 
In Fig. |6(b)] we have shown two different k < solutions, one for v = —2.0 and the other one for v = 2.0. Here, we 
see that, for both solutions, the Ricci scalar is symmetric, negative everywhere and asymptotically AdS. Additionally, 
it can also be shown that, it is asymptotically flat for large negative value oi v {v <^ —1) or very small positive value 
of (0 < ^ 1). In these solutions we also note the energy density to be symmetric as I — > — L For v = 2.0, though 
Pt is positive for all it has a minimum value at / = and it increases for |/| > tending to a constant value for 
large 

IV. CONCLUSIONS 

The main purpose behind this paper was to find simple solutions of EiBI gravity in three dimensions. To this end 
we have found an analytical FRW cosmology for pressureless dust and analytical circularly symmetric solutions. We 
have also obtained the numerical solution for radiation, in the case of cosmology. Apart from a flat fc = spatial 
section, we have also explored the k — +1 and A: = — 1 FRW models. The cosmological solutions presented here seem 
to point to a generic feature that singularities are not present if k < whereas singularities do arise if k > 0. It is 
not clear whether one can make any general statement on this based on the solutions we have found. In the class of 
circularly symmetric, static solutions we note that they are nonsingular for both k > and k < 0. 
One may also look for more general solutions in the circularly symmetric case, with a non-constant goo- This is likely 
to lead to three dimensional black holes in EiBI theory. Finally, it will surely be useful to see if our solutions in 
this toy three dimensional theory can serve as viable pointers towards finding new analytical or numerically obtained 
geometries in the real four dimensional scenario. 
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FIG. 6: Comparison of solutions with different energy-density and curvature scalar, for various values of the parameters (/x and 
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